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Abstract 

The general relativistic effects on spin-flavor oscillations above the core of 

type II supernovae are investigated. The evolution equation is derived and 

the relative magnitudes of the terms in the Hamiltonian, which arise from 

the weak, electromagnetic and gravitational interaction, are compared. The 

effects on the resonance position and the adiabaticity are studied. Explicit 

calculations are presented for non-rotating and slowly rotating stars. 
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I. INTRODUCTION 



If the neutrino possesses a magnetic moment, it has been shown that its spin can be 
flipped while passing through a magnetic field [|l|-[5|]. This can lead to spin-precession and 
spin-flavor oscillations. In the standard model of particle physics the neutrino magnetic 
moment is zero. However non-standard neutrino properties such as mass and magnetic 
moment attract enormous attention since they are not only natural by-products of many 
grand unified theories; but they can also lead to neutrino oscillations, which provide a very 
popular and elegant solution to the long-standing solar neutrino problem M. 

Terrestrial experiments constrain the magnetic moment of the neutrino to values < 
lO -10 ^ QHH? where /ie is the Bohr magneton. Astrophysical bounds are more stringent 
constraining fi < 2 x 10 _12 /i^ jToHl^ll- The latter are, however, somewhat model-dependent. 

Thus for spin-flips to be discernible, very strong magnetic fields are required. One of the 
places where magnetic fields of such strengths are likely to exist is above the proto-neutron 
star of a type II supernova (SN). Observations of magnetic fields on the surface of white 
dwarfs suggest the existence of magnetic fields of up to 10 10 G at the surface of the iron core 
15| ; energetic considerations suggest even higher fields [|16j. 

In type II SNe neutrinos play a crucial role since they are believed to drive the supernova 
explosion |T^,|T3. Neutrino oscillations could therefore affect the explosion mechanism di- 
rectly. Also neutrinos from a nearby SN explosion, SN 1987 A, have been detected on Earth 
[IJ§ , and the measured neutrino spectrum has been used to test various oscillation scenarios 



UH. Neutrino spin-flips have also been invoked to explain the high pulsar kick velocities, 



as the interaction with a magnetic field can lead to an asymmetric neutrino emission p2 |. 



It is known for some time that gyroscopes precess in strong gravitational fields [23,24 



Here I will investigate under which conditions gravitational effects can influence the neutrino 
spin flip in a SN. 

SN explosions occur at the end of the life of massive stars and leave a neutron star 
or a black hole behind. Therefore it is of importance to study neutrino oscillations in 
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curved space-time and calculate the magnitude of general relativistic effects. Cardall & 
Fuller have calculated the gravitational effect on neutrinos in vacuum and matter. 
They derived MSW (Mikheyev-Smirnov-Wolfenstein) oscillation formulae in curved space- 
time and presented expressions for radial and azimuthal propagation in a Schwarzschild 
geometry. Piriz et al. [2B| derived the resonance conditions for neutrino conversions with 



regard to neutrino oscillations in the vicinity of active galactic nuclei. They present explicit 
calculations for arbitrary orbits in a Kerr-geometry. 

Here I will briefly rederive the evolution equation for massive Dirac neutrinos interacting 
with a magnetic field in matter. In Sec. 3 and 4 I will study how the neutrino conversion 
is affected by the curvature of space-time near a static and a slowly rotating proto-neutron 
star. 

Throughout this paper I shall set G = h = c = 1. 

II. THE EVOLUTION EQUATION IN CURVED SPACE-TIME 

Gravitational effects on the spin arise through the "spin connection" appearing in the 
Dirac equation. This effect has been studied by Cardall & Fuller [25] and somewhat more 



formally by Piriz et al. |26| (also see [27,28|). In the following I will follow the treatment 



of Ref. ||25|| . The Dirac equation in curved space-time (ignoring magnetic fields and matter 
effects for the moment, and restricting myself to two neutrino species) for massive neutrinos 



interacting minimally with the gravitational field can be written as follows |29] 



[ 7 °eS(a /l + r |l ) + M]v = o. (l) 

Here ip is a column vector of spinors of different mass eigenstates, the 7 a 's represent the 
Dirac matrices and M is the mass matrix, M 2 = diag[mf,m|], where mi and m-i are the 
mass eigenvalues of the two neutrino species. I follow the convention that greek indices 
refer to general curvilinear coordinates, while the latin indices a, b, c, d refer to locally in- 
ertial (Minkowski) coordinates. The tetrads connect these sets of coordinates. g^ v will 
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denote the metric of curved space-time and rj a b the metric in flat space, where I will use the 
convention that i] ab = diag[— 1, 1, 1, 1]. The tetrads then satisfy 

9»v = Vabefel (2) 

The explicit expression for T M is 

T, = \[l\l c Ke cu ,,. (3) 

Here the semicolon denotes a covariant derivative and a comma an ordinary derivative. It 
can be shown that 

71 = 2r/V " 2rT 7 b - 1ie dahc lbld) (4) 

where e abcd is the (flat space) totally antisymmetric tensor, with e 0123 = +1. With Eq. (|]), 
the non- vanishing contribution from the spin connection is 



7 a e£r M = 7 X [iA GfM 



75 

" 2 



(5) 



where 



A% = \e»e abcd (e bl/ , a - e b ^)e^ d . (6) 

The expression in Eq. (|5|) treats left- and right-handed states differently. In order to group 
it with terms arising from matter effects, we can without physical consequence add a term 
proportional to the identity to obtain 

7^ = 7^^^), (7) 

where Vl is the left-handed projection operator Vl = (1 — 7 5 )/2. As pointed out by Cardall 



& Fuller P5 |, is diagonal in spin space and cannot produce spin flips on its own. However 
it can change the resonance as well as the adiabaticity conditions in the presence of other 
off-diagonal terms such as those arising from an interaction with a magnetic field. 

A further complication arises in curved space-time: The neutrino trajectories are no 
longer necessarily straight lines as in flat space-time. They can be written as curves of the 



affine parameter A. The tangent vectors to the null worldline, p M = dx^/dX, can be found 
by solving the geodesic or Hamilton- Jacobi equation (see, e.g., P(|). As suggested in Ref. 
|25| the four-momentum is then constructed by assuming that the neutrinos are energy 
eigenstates and that the three momentum is collinear with {pi,P2,P3)- 



The mass matrix M is related to the mass matrix in flavor space Mf through a rotation 
by the vacuum mixing angle 9 V as follows: 

(m\ \ 

M} = U[ )U\ (8) 

V mi) 

with 

(cos 9 V sin 9 V \ 
(9) 
— sin 9 V cos 9 V ) 

The neutrinos interact with matter through a weak interaction current of the form: 
A^Vl, where in the chiral basis {veL^eR^^^R)- = diag[V c + V n , V n , 0, 0] (Here the z^s 
can, of course, be replaced by z/ r s.). The matter potentials due to the neutral and charged 
current in neutral, unpolarized matter are given by V n = —\[2GFn n and V c = 2^/2Gpn e , 
respectively, where Gp is Fermi's constant and Ei the energy as measured in a local inertial 
frame. n e and n n are the electron and neutron density in the rest frame of the background 
matter. 



The interaction of the neutrino with an external magnetic field is described by a term of 
the form 

fiv ab F ab ^ } (10) 
where fi is the magnetic moment matrix, F a b the electromagnetic field tensor and a ab = 

ih a ,i b }- 

Only considering conversions between two neutrino flavors we can write the evolution 
equation as follows: 
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(11) 



where H is the Hamiltonian in the chiral basis. Including the weak, electromagnetic and 
gravitational interaction the Hamiltonian can be written as: 
/ 



H = 



V c + V n + P IX A G /E l - cos 29 v 



sin 2^ 

HeeB 



fi ee B 

V n + P iX A G jE l + cos 29 v ^ e B 



\l>eiiB 
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-gcos2# ? , ^#sin2#„ 
fgsin2^ ffcos20„ 



„2 



(12) 



where B is the component of the magnetic field which is perpendicular to the neutrino 
trajectory and Sm 2 the difference of the squares of the neutrino mass eigenvalues in matter: 
Sm 2 = m\ — m\. Here it was assumed that the right-handed components are sterile and do 
not interact with matter. Resonance occurs whenever the diagonal elements cross. From 
the expression for H we can immediately read off the resonance conditions for spin-flips and 
MSW-conversions: 
For v eL <-> v^: 



for v eL <-> v eR : 



5ui 2 

V c --—coa26 v = 0, 



V c + V n + P ll A G jE l = 0, 



(13) 



(14) 



for U eL <r+ V^R. 



V c + V n + P IX A G /E l 



Sm 2 

m 



cos 29 v = 0, 



(15) 



for is^l ^ v eR : 



V n + P^/Et + — r cos 26 v = 0, (16) 
ihji 

and for v^ L <-> z/ Mfl : 

K + = 0. (17) 



These equations are similar to those derived in Ref . [ 26j . In the active-active MSW transition, 



Eq. (|13D, the matter potential is due to the charged current only and proportional to Y e . 
[Y e is the electron fraction: Y e = n e /(n n + n e ).] However for the transitions corresponding 
to Eqs. ( p!4|) and (p~5f ) the potential is due to the neutral as well as the charged current and 
thus proportional to (1 — 2Y e ). The latter potential is strongly suppressed in the isotopically 
neutral region (number of protons=number of neutrons) |33|| . 

In the following sections I will calculate the 'gravitational current', P^Aq, for a static as 
well as for a rotating star. 



III. RELATIVISTIC EFFECTS NEAR STATIC STAR 



The geometry in a spherically symmetric, static spacetime can be globally represented by 
the Schwarzschild coordinate system =^ (t, r, 6, 0). MSW transitions in a Schwarzschild 
metric have already been discussed in Ref. |25] and their findings can easily be extended to 
spin-flip conversions. The Schwarzschild line element, which serves to define these coordi- 
nates, is 



ds 2 = g^dx^dx" 



_ e 2*(r)^2 + e 2A(r) rfr 2 

+r 2 d0 2 + r 2 sin 2 ■ 



(18) 



Using the tetrads 



diag 



e -*(r) e -A(r) 



1 



r r sin 9 



(19) 
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it is easy to verify that Aq = 0, as pointed out by p5 |. In Schwarzschild metric it is common 



to write e 2 ®^ = l — r s /r and e 2A( - r * ) = 1/(1 — r s /r), where the Schwarzschild radius r s = 2A4 
and M. is the central gravitating mass. The energy-like component of the four-momentum 
is the energy and a conserved quantity: Pq = —E. The energy measured by a locally inertial 
observer is Ei = Ee^®^] i.e., the energy of the neutrino suffers a gravitational redshift by 



l/Jl — r s /r [p4 |. Thus the resonance conditions are merely altered by setting E\ equal to 



Ee~^ r \ and Aq = 0. 

For a particular transition one can define a helicity mixing angle, 6 m , which diagonalises 
the corresponding submatrix in the Hamiltonian (see, e.g., ||). For instance, for the tran- 
sition v eL <-> v^r, the mixing angle is defined by: 

tan 29 m = ^ _ ^ - . (20) 

The adiabaticity parameter corresponding to a particular transition, 7, quantifies the mag- 
nitude of the off-diagonal elements compared to the diagonal ones of H in the basis of the 
instantaneous eigenstates. At the resonance point it must be greater than 1 for adiabatic 
propagation. 7 is modified by relativistic effects, and for the case v eL <-> u^ R) can be written 
as: 

, x 8e-*^Efj? eu B 2 , 

ry( r ) = Lit 21) 

n res) 5m 2 cos26 v e- A ^)\d\n(V c + V n )/dry 1 ; 

where the logarithmic derivative is evaluated at the resonance position. Thus the adiabaticity 
of the evolving neutrino is affected by the spatial dependence of the metric, which appears 
for example in the determination of the local energy. The conversion probability can then 
be expressed as 

Pu eL ^ R = \-(l- P ho P ) cos 26 v cos 26 m , (22) 

where the non-adiabatic hopping probability is given by Ph op = exp — |7(r res ) ■ 

The conversion probabilities corresponding to the other transitions of Eqs. ( P^ - |T7| ) can 
be found accordingly. 
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IV. RELATIVISTIC EFFECTS NEAR ROTATING STAR 



Rotation of the collapsing star will cause a dragging of the inertial frames and will result 
in off-diagonal terms in the metric g^ u . These off-diagonal terms can now flip the spin of 
the neutrinos by transferring the angular momentum of the star to the neutrino spin. Using 
perturbation theory one can write down a metric to first order in the angular velocity of 
rotation fl of the central mass (e.g., see ([H],[J^]): 

ds 2 = _ e 2*(r) dt 2 + e 2A(r) dr 2 

+r 2 (d6 2 + sin 2 Odcf) 2 ) - 2ur 2 sin 2 Od^dt, (23) 
where u is the local precession frequency and can be written as u = 2J/r 3 [31], where J is 



identified with the angular momentum of the central gravitating mass. Here it is assumed 
that the angular momentum vector of the star, J, is aligned with the polar axis. The metric 



of Eq. (|23|) is valid as long as the rotation is slow, i.e., for ila « 1, where a is the radius 
of the central mass. 

In this metric purely radial motion is only possible along the axis of rotation as expected 
from symmetry considerations. From these same considerations it is also evident that motion 
in a plane is only possible in the equatorial plane. 

However, here I will confine myself to nearly radial trajectories and work in the weak- 
field limit so that the curvature of the neutrino path can be neglected. Thus one finds from 



Eqs. (g) and (H): 

P„AyEi « e~ A(r) Jr ~ 3 cos 9, (24) 

for relativistic neutrinos moving radially outwards (i.e. P\ ~ Pq). It can be noted that 
for propagation in the equatorial plane the 'gravitational current' is still zero. But at any 
other polar angle there will now be an additional term in the Hamiltonian which will shift 
the resonance position and modify the adiabaticity conditions. In Eqs. (|20| ) and fl2~T|) one 
can then simply replace V c + V n by V c + V n + P^Aq/Ei. It is interesting to note that P^Aq 



changes sign at the equator and could thus lead to an asymmetry in the neutrino emission. 
It is evident from Eq. fl2^]) that the 'gravitational current' moves the resonance point closer 
to the center of the SN in the upper hemisphere and farther away in the lower hemisphere. 
The magnitude of this effect will be discussed in the following section. 



V. DISCUSSION 

Simple angular momentum arguments suggest high angular velocities for neutron stars, 
which have also been verified by observations. During collapse to a neutron star the angular 
velocity of a typical star will increase from about Qo ~ 10 _6 s~ 1 to Q ~ 10 4 s _1 . Fig. 1 shows 
the absolute magnitudes of the various components in the Hamiltonian given by Eq. ( |i"2|) for 
a neutrino propagating in a typical SN environment. I consider here the case when the proto- 
neutron star is rotating, as Aq = otherwise. The neutrino is assumed to move radially 
along the polar axis. The matter potential was calculated using the progenitor model for 



a 15 M Q , low metallicity star as calculated in Ref. [O] on the basis of the models of Ref. 



35| . It can be seen that in the region between 0.02 R Q and 0.3 R@ the matter potential is 
strongly suppressed due to the high abundance of isotopic elements such as Oxygen, Carbon 
and Helium. 

The magnetic field was assumed to be a dipole field which is aligned with the rotation 
axis and is given by: 

B = B Q (— V sintf, (25) 



r 



where Bq is the magnitude of the magnetic field at the surface of the iron core at a radius 
r$. The value of the neutrino magnetic moment was taken to be /i = lO" 12 ^ and -Bo was 
set to 10 10 G. The 'gravitational current' was calculated assuming that the iron core has 
collapsed to form a proto-neutron star of mass M = 2M Q and a radius of a = 10 km. Its 
angular velocity was taken to be = 10 V 1 . It is seen that the 'gravitational current' 
comes relatively close to the matter potential in the isotopically neutral region of the SN 
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envelope. Also shown in Fig. 1 is the absolute difference of the diagonal elements of the 
Hamiltonian for the transition v e i <-> v^r, i.e. the right-hand side of Eq. fll4l). Resonance 
occurs whenever the difference in the diagonal elements goes through zero. It is evident that 
in order for the 'gravitational current' to affect the neutrino spin-flip, the resonant transition 
has to take place in the isotopically neutral region. For a choice of 5m? = 10 -10 eV 2 and 
cos^ ~ 1 (large mixing angle) the resonance occurs at about r = 0.15i? Q . 

For magnetic fields as high as indicated in Fig. 1 the off-diagonal term, fiB, dominates 
the potential by several orders of magnitude such that in the isotopically neutral region 
spin-flavor precession of almost maximum depth occurs. Strong precession occurs whenever 
\iB is much greater than AH and the precession amplitude is given by ||16|| : 

A - ^ B)2 (26) 

In that case the gravitational current has a negligible effect on the conversion probability. 
However for weaker - and possibly more realistic - magnetic fields (B < 10 7 G) [iB becomes 
comparable to AH, suppressing the spin-flavor precession. However [iB is still greater than 
AH at the resonance layer, where resonant conversion occurs as long as the adiabaticity 
condition is satisfied. 

It is readily seen that for our example the 'gravitational current' is still by about one 
order of magnitude smaller than the matter potential and will thus not significantly affect 
the resonance position or the adiabaticity in the envelope of the SN (even for B < 10 7 G) . 

However if the mass of the proto-neutron star or its rotation rate is increased, or the 
matter potential in the envelope is suppressed by about one order of magnitude (i.e., by 
reducing the abundance of high Z metals), the presence of the gravitational current will 
start to shift the resonance position and affect the transition probability noticably. 

The conditions for which the gravitational contribution becomes important are neither 
contrived nor unrealistic. These calculations are only intended to yield crude orders of mag- 
nitude estimates for the different terms in the evolution equation and more exact numbers 
would require a more careful modelling of the collapsing star. For instance, I have neglected 
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the infall velocity of the matter during the collapse and only considered radial neutrino 
trajectories in the weak- field limit. However these simple estimates presented here show 
that the terms arising from the space-time curvature can potentially become significant for 
spin-flip transitions in the SN envelope even though they are probably not important in 
most realistic cases. The extent to which they affect the neutrino transition depends on the 
detailed composition of the SN, as well as the mass and rotation rate of the proto-neutron 
star. 
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FIGURES 



FIG. 1. Magnitudes of the various components in the Hamiltonian of (12) as a function of 
distance from the center of a 15 Mq type II SN. The energies are measured in eV and the radius is 
measured in units of solar radii. The lines represent the matter potential (solid), the gravitational 
term (short-dash) , the magnetic term (dotted) and the absolute difference of the diagonal elements 
in the Hamiltonian (long-dash). 
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